Algebraic and Analytic Properties of the 
One-Dimensional Hubbard Model 

Frank Gohmann^ and Shuichi Murakami-^| 

^Department of Physics, Faculty of Science, 

University of Tokyo J]] 
Hongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 
■f Department of Applied Physics, Faculty of Engineering, 
University of Tokyo, 
Hongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 

Abstract 

We reconsider the quantum inverse scattering approach to the one-dimensio- 
nal Hubbard model and work out some of its basic features so far omitted in 
the literature. It is our aim to show that i?-matrix and monodromy matrix of 
the Hubbard model, which are known since ten years now, have good elemen- 
tary properties. We provide a meromorphic parametrization of the transfer 
matrix in terms of elliptic functions. We identify the momentum operator 
for lattice fermions in the expansion of the transfer matrix with respect to 
the spectral parameter and thereby show the locality and translational in- 
variance of all higher conserved quantities. We work out the transformation 
properties of the monodromy matrix under the su(2) Lie algebra of rotations 
and under the //-pairing su(2) Lie algebra. Our results imply su(2)©su(2) in- 
variance of the transfer matrix for the model on a chain with an even number 
of sites. 
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1 Introduction 

The one-dimensional Hubbard model is one of the most thoroughly studied integrable 
quantum systems with applications in solid state physics. Starting with the seminal 
article of Lieb and Wu lots of its physical properties have been worked out exactly 
P]. For the case of half-filled band, in particular, a complete picture of its elementary 
excitations is available by now 0, All excited states are scattering states of only 
four quasiparticles, two of which carry spin but no charge, whereas the other two carry 
charge but no spin. The S'-matrix of these quasiparticles has been calculated exactly. 
These achievements give a precise meaning to the notion of spin-charge separation in 
one-dimensional solids. 

All exact results on physical properties of the one-dimensional Hubbard model ob- 
tained so far rely on the extensive use of the coordinate Bethe Ansatz. Since Bethe 
wave functions, however, are difficult to handle, any exact calculation of local quantities 
going beyond the long distance asymptotics of correlation functions 0] seems to demand 
for an algebraic treatment. Moreover, an algebraic treatment is likely to facilitate the 
calculation of the thermodynamical properties of the model 0. 

At present we know two algebraic structures related to the Hubbard model, a graded 
Yang-Baxter algebra, developed in the works of Shastry |^, ^, Q and Olmedilla et al. 
TT| , p!2| and a representation of the Y(su(2)) Yangian quantum group commuting 



with the Hubbard Hamiltonian, which was discovered by Uglov and Korepin |TB[. The 
relation of these two notions was recently exposed by the authors |T^. 

Although i?-matrix and L-matrix of the Hubbard model are known since long, it 
took nearly ten years before it was shown that the i?-matrix satisfies the Yang-Baxter 
equation An algebraic Bethe Ansatz was performed only recently in a remarkable 



preprint by Ramos and Martins [|T^ . Every progress in the development of an algebraic 
approach was hindered before by the complexity of i?-matrix and monodromy matrix 
and by several unusual features of these basic tools of the quantum inverse scattering 
method. The monodromy matrix is 4 x 4 rather than 3 x 3, as one might have guessed 
naively from the fact that there are two levels of Bethe Ansatz equations. It further 
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seems to be impossible to find a parametrization of tlie i?-matrix, sucli tliat it becomes 
a function of the difference of tlie spectral parameters. 

Amazingly, not even the most elementary properties of i?-matrix and monodromy 
matrix have been worked out so far. For that reason we look again at the construction of 
the graded Yang-Baxter algebra. We begin with a description of the spin model [0) 01 
that is related to the Hubbard model by means of a Jordan- Wigner transformation in 
section 2. Our account is based on the Yang-Baxter equation. In section 3 we show that 
there exists a meromorphic parametrization of the transfer matrix in terms of elliptic 
functions. Section 4 is devoted to the rederivation of the graded Yang-Baxter algebra 



11 . We show how to treat general twisted boundary conditions. As a byproduct 



we find a simple method to obtain higher conserved quantities of the Hubbard model 
from their counterparts for the spin model. These conserved quantities are generated 
by the graded trace of the fermionic monodromy matrix. We identify the momentum 
operator for fermions in the zeroth order of the expansion of this generating function with 
respect to the spectral parameter. Thus all higher conserved quantities are local and 
translational invariant. In section 5 we derive the properties of the monodromy matrix 
under a combined particle-hole and gauge transformation, which is characteristical for 
the Hubbard Hamiltonian. It turns out that the graded trace of the monodromy matrix 
is invariant under this transformation, if the model is considered on a chain with an even 
number of sites. In section 6 we investigate the behaviour of the monodromy matrix 
under su(2) transformations. Our results are complementary to the work of Ramos 
and Martins [|T^ and should provide a means to discuss the symmetry properties of 
quasiparticles within the algebraic approach. Like the Hamiltonian, the graded trace 
of the monodromy matrix turns out to be invariant under rotations of the spins and. 



if we consider an even number of sites, also under the r^-pairing su(2) Lie algebra fl^, 



18| , [19| , PQ] . In order to make our representation self-contained, we expose the i?-matrix 



along with a list of relations among its elements in appendix A. Appendix B provides a 
detailed discussion of the momentum operator on a lattice. In appendix C it is shown 
how to obtain the momentum operator from a monodromy matrix for free fermions. 
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2 The spin model 

As all one- dimensional fermionic models the Hubbard model is related to a certain spin 
chain by a Jordan- Wigner transformation. The Yang-Baxter algebra corresponding to 
this spin chain is easier to formulate and closer to intuition than the graded Yang-Baxter 
algebra of the Hubbard model. In fact, Shastry in his seminal articles H was using 
the language of spins rather than the language of electrons. The graded form of the 
Yang-Baxter algebra having the advantage of being defined directly in terms of fermi 



operators was derived later by Olmedilla et al. [|T^, [TT], |12|. Let us follow the historical 



route here and start with a description of the spin model. Its Hamiltonian is 

L 
J = l 

Here and in the following we are using implicit summation over doubly occuring indices. 
H describes a periodic spin chain of L sites (cr^_|_i^ := crf^) with two species of spins, 
labeled t and |, at each site. U is the strength of an on-site Ising coupling between 
the species. The interaction between nearest neighbours is of XX-type for both species 
independently. Thus, in the limit t/ — the model decouples into a pair of non- 
interacting XX-chains. 

As was shown by Shastry, H is the logarithmic derivative of a certain transfer matrix, 
which can be obtained by appropriately coupling together two copies of the Yang-Baxter 
algebra of the XX-chain. The i?-matrix of the XX-chain is 

r = i(a + 6+(a- 6)(t^ ® a' + a"" a"" + ® a^) , (2) 

where a and b have to satisfy the free fermion condition 

a^ + b^ = l . (3) 

If we introduce the parametrization a = cos(A), b = sin(A), then r = r(A) satisfies the 
Yang-Baxter equation in the form 



ri2(A - /u)ri3(A)r23(Ai) = r23(Ai)ri3(A)ri2(A - /i) (4) 
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which is sometimes called difference form of the Yang-Baxter equation. r(A) is regular, 
i.e. r(0) = P, the permutation of the two factors of C^. As usual, the indices in 
refer to the canonical embeddings of r(A) into (S> C^. The L-matrix /(A) of the 
XX-chain is the fundamental representation /(A) = r(A) of the Yang-Baxter algebra, 

f(A-/x)(/(A)®/(/x)) = (/(/x)®/(A))f(A-/x) , (5) 

generated by f(A) := Pr(A). As usual, /(A) in eq. (j^) has to be understood as matrix in 
auxiliary space with entries acting on a quantum space. 

For the construction of R- and L-matrices corresponding to the spin Hamiltonian (|l]) 
we have to duplicate the above construction by attaching a label referring to the spin 
species to each cx-matrix. The matrix r(A) is replaced by r,-(A) (r =|, |), and we may 
redefine r(A) as 

r(A):=r^(AK(A) . (6) 

This new i?-matrix obviously satisfies (§) and is again regular in the sense that r(0) = 
P^Pj =: P is a permutation operator. Whenever an explicit matrix representation is 
required we will use the convention a" = cr" ® /2, cr° = J2 ® a", where I2 denotes the 
2x2 unit matrix. Correspondingly, the 4x4 unit matrix, which will be needed below, 
is denoted by 14. Now Shastry's P-matrix associated to the Hamiltonian reads 

P(A,/x|/i,/) := ch(/i-0^%-^+ sh(/i-/)^%ti^a^®a^®/4 . (7) 
^ ' ^ ^ ^cos(A-/i) ^ ^cos(A-F/x) ^ ' 

This is a four parameter family of 16 x 16-matrices. It was shown by Shiroishi and 
Wadati |T^], that it satisfies the Yang-Baxter equation in the form 



Pi2(A, yul/i, /)Pi3(A, v\h, m)P23(/i, z/|/, m) = P23(yU, m)Pi3(A, v\h, m)Pi2(A, /) , 

(8) 

provided that the parameters are constrained by the equations 

U sh(2/i) sh(2/) sh(2m) 



4 sin(2A) sin(2/i) sin(2z/) 
This constraint can be satisfied for arbitrarily small m and v. Hence, 



(9) 



Ljk{\\h) := cos(A)Pjfc(A, 0\h, 0) = r,fc(A)e'^'^lTSi (10) 
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is a representation of the Yang-Baxter algebra generated by the i?-matrix (|^. The 
index j in (|T0|) refers to the auxihary space, the index k to the quantum space. Note 
that Ljk{X\h) 7^ Lkj{X\h). Ljk{X\h) is a tensor product of two XX-chain L-matrices 
coupled in auxiliary space. If we solve the first equation in (Q) for h and choose the 
branch of solution properly, then h{X = 0) = 0, and Ljk according to (p!0| ) is again 
regular. For U = the constraint (|^) is satisfied by /i = for all A. We get back the 
free model as it should be by construction, and the i?-matrix (0) becomes a solution of 
the Yang-Baxter equation in difference form (Q). On the other extreme, we can carry 
out the limit f/ — oo after properly rescaling A X/U. Then the constraint @ 
becomes A = 2sh(2/i). However, A disappears from the definitions of Lj^ and Rjk, since 
cos(A/t/) — > 1 and sin(A/f/) — > for every A. Thus rjk{{X — jj)/U) Pjk, and 

i?,fc(A,/i|/i,/)^L,fc(0|/i-/)=P,,e("-')'^lTSu . (11) 

It is easily verified that the expression on the right hand side of ([TT|) satisfies the Yang- 
Baxter equation in difference form (^) and that the corresponding Hamiltonian is the 
on-site part of (|l]) with U = 1. Hence A is the natural spectral parameter of the model at 
zero coupling, whereas h is the natural spectral parameter in the strong coupling limit. 
The algebra (|^), (|), (|]) interpolates between these limits. For the remainder of this 
article we will suppress the arguments h and / of the R- and L-matrices, assuming that 
they are given as functions of A and ^ by the constraint (P). 

To finish the description of the spin model let us introduce its monodromy matrix, 

Tl(A) :=L,z.(A)...L,i(A) . (12) 

The index a in this definition refers to the auxiliary space. The transfer matrix of the 
spin model, t{X), is the trace of Tl{X) over the auxiliary space. It follows from the 
regularity of Ljk that U := t{0) is the shift operator for spins, 

Uaf^ = af^,,U , J = 1,...,L , r =t, i . (13) 

A brief calculation yields the derivative of the L-matrix at zero spectral parameter. 
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This equation implies that the Hamiltonian is obtained from the transfer matrix t{X) 
as logarithmic derivative, H = d\ ln(t(A))|^^g. 



3 A meromorphic parametrization of the transfer 
matrix 

The considerations in this section were motivated by two facts. First, in case of the eight 
vertex model Baxter's meromorphic parametrization of the i?-matrix yields a solution 



of difference form of the Yang-Baxter equation [^. Second, the transfer matrix en- 
ters certain functional equations, the solutions of which usually require strong analytic 
properties. 

Let us solve the constraint (|D for e"^^, 

2^ [/sin(2A) r f/2 sm\2X) 



(15) 



4 V 16 

The only possibility to remove the square root on the right hand side is by replacing 2A 
by am(M), the amplitude function and setting k = ill /A. Then e^^ is expressed in terms 
of Jacobi elliptic functions as 



^2h 



dn{u) — ik sn{u) 



(16) 



Because of the homogeneity of the Yang-Baxter algebra, we may multiply Ljk{u) by e^. 
Then h{u) enters into the definition of the monodromy matrix Ta{u) only as meromorphic 
function e^^ of the redefined spectral parameter u. Letting A := {a + 6)/2 and B = 
{a — b)/2 we see that the matrix rjk in the definition of the L-matrix is of the form 



^ e o o e ^ 



o e e o 
o e e o 
e o o e 



(17) 



where e denotes an even polynomial in A, B and o denotes an odd one. e and o 
are operators on quantum space k, whose precise form is irrelevant for the following 
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arguments. The rules = = e, eo = oe = o, e + e = e and + = imply 



^ e o o 



o e e o 
o e e o 
e o o e 



^ e o o 



o e e o 
o e e o 
e o o e 



(18) 



This means that the monodromy matrix (|T2D is again of the form (0). Therefore t{u) 
is an even polynomial in A and B. In other words, t{u) is a polynomial in A^, and 
AB. Since 



A' 



1 + sn(u))/4 , 5^ = (1 - sn(M))/4 , AB = cn(u)/4 



(19) 



t{u) is a meromorphic function of u. 

To state the problem of the analytic structure of the i?-matrix let us go one step 
back and write again a for cos(A) and b for sin(A). Let furthermore c := e^^. Then a, 
b and c are connected by the free fermion condition (|]) and the constraint (|^), i.e. a, b 
and c are lying on a complex curve given by the algebraic equations 



+ b^ = 1 , c - 1/c = Uab 



(20) 



This curve may be called the spectral curve of the Hubbard model. Unfortunately we 
could neither find a meromorphic parametrization of this seemingly simple structure nor 
assign a geometrical meaning to it. 



4 The Hubbard model 

Applying a Jordan- Wigner transformation to R- and L-matrix of the spin model we 
obtain the graded Yang-Baxter algebra |22] of the Hubbard model [|lT|. No extra ef- 



fort is necessary to introduce the grading. It is rather induced by the Jordan- Wigner 
transformation. Before turning to the fermionic formulation let us perform a gauge 
transformation of R- and L-matrices with 4x4 transformation matrix 

^^_gX(aW)/2 _ ^21) 
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Then 

L,(A) LkiX) := GMX)G^' = Ghilk^{X)(^ki{X))Gh , (22) 

i?(A,/i) i?(A,/i) := {Gh(E)Gi)R{X,^i)iG^^(g)Gr^) . (23) 

We suppress the auxihary space index of the L-matrix here and in the following and 
consider I/fc(A) as 4 x 4 matrix with entries acting on quantum space k. Unlike -R(A, fi), 
the transformed matrix -R(A, /x) is symmetric. By use of the definition 

R{X,ij) := PR{X,i2) (24) 

the Yang-Baxter algebra assumes the form 

i?(A, /x) (l,(A) ® L,(/i)) = (L,(/i) ® Lfc(A)) RiX, fi) , (25) 



which is most convenient for changing to a fermionic formulation 

Since the canonical anticommutation relations for fermi operators are invariant under 
gauge transformations and under particle-hole transformations, there is some freedom 
in the definition of the Jordan- Wigner transformation. Using the abbreviations 

k k 

Uk-=J2^n ' (^k-=J2^jl ,k = l,...,L, (26) 

i=i 3=1 



where the electron densities and setting further uq := do := 1, we take the choice 

'^t^ = 4/'''''-' , a,-^ = 0,^6-^'^'^'=- , (27) 
ati = 4/'''^'''+''-'^ , a,-^ = Cfcie-''^("-+'^'=-) . (28) 
The XX L-matrices /fc|(A) and lki{X) can now be expressed in terms of fermi operators, 

4-^(A) = e-^""'="'/2£^^(^)gi7r«,_iaV2 ^ (29) 

4^ (A) = e--(«i+d.)-V2£^^(^)gi-(«i+rf.-iV/2 ^ (3Q) 
where CkriX) (r =t, |) is defined as 

, Sin(A) + id^Hkr Ckr 

A.(A):=| I . (31) 

-icj^ cos(A) - e'^Ukr 
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£fcT-(A) is an L-matrix for free fermions. Its entries for different quantum space indices 
k either commute or anticommute. This fact can be formally described by assigning 
a parity TrlCkri^Yj) = 0, 1 to each matrix element. Call a matrix element even, if its 
parity is zero, odd if it is one. Odd matrix elements with different quantum space indices 
anticommute, whereas even elements commute with all elements with different quantum 
space indices. A grading is a function p, which assigns parity to the basis vectors of the 
auxiliary space. Let p{i) := p(ej). Then, with the grading p{l) = 0, p(2) = 1, p E Z2, 
the elements of the free fermion L-matrix have parity 7r(£fc^(A)*) = p{i) +p{j). The 
crucial point about these notions is that they allow for the introduction of a graded 
tensor product ® which respects comultiplication. Let ® be defined by the equation 

(Af e)^' := . (32) 

Then f f = (AC f BD) for all matrices A, B, C, D of the parity defined 

above. The graded tensor product can be used to formulate a graded Yang-Baxter 



algebra p2[. Grading and graded tensor product may be defined for arbitrary matrix 
dimensions. 

In appendix C we give a brief account of the free fermion model which will be needed 
below to recover the lattice momentum operator from the monodromy matrix of the 
Hubbard model. The free fermion monodromy matrix is 

TLri\):=CLA\)...Cir{\) , (t =T, i) • (33) 

Note that our choice of Jordan- Wigner transformation and thus our free fermion L- 
matrix differ from that in |lT[ . 

Inserting (|2S|), (|§) into (0) we find 

Lk{X) = WV,''Ck{X)Vk-iW~' , (34) 



where 



W := diag(l,l,i,i) , (36) 
A(A) = G,,{Ck^i\)®Ckii\))Gh . (37) 
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The grading comes in naturally in (p?]), since the operator I2 ® g-i'r(«L+'ife)o-^/2 ^\^qt^ 
moved to the left in the tensor product /fct(A) (8'/fcj(A) induces a gauge transformation on 
£fcl(A), which affects only the odd elements. We will see below that Ck{\) according to 
equation (^) is an L-matrix for the Hubbard model. The associated /2-matrix follows 
from (|25|). First of all we have 

Lk{\) ® LkijJ,) = 

{W ® W) (V,-' ® VjT') X (Ck{X) f X-^ (\4„i ® [W"' ® W~^) , (38) 

where ® is a graded tensor product of 4 x 4 matrices (cf. (^)) with grading p{l) = 
p(4) = 0, p{2) = p{3) = 1. We are using the same symbol for graded tensor products of 
2x2 and 4x4 matrices. For 4x4 matrices the grading will always be as above, and 



for 2x2 matrices we will use the grading introduced below (32). The matrix X in (pSf) 
is the diagonal matrix 

X:=a"®diag(l,i,i, 1)® J2 . (39) 



Let 



r := diagfe 



ia \p i7 A5 



e") , 



where a, j3, 7, 5 may generally be mutually commuting operators. Then 



^(A,/i),r® F 



^ q; + 5 = /5 + 7 mod27r 



Since Vk and W are of the form (]40|) and satisfy (|4^) , we may infer from (^) that 



7^(A, /i) [Ck{\) f A(/i)) = (A(/x) f A.(A)) 7^(A, /i) 



(40) 



(41) 



(42) 



(43) 



where 



7^(A,/i) = X-i^(A,/i)X . (44) 
Hence, the L-matrix of the Hubbard model is a representation of the graded Yang-Baxter 



algebra. This result is due to Olmedilla et al. [jl^]. To be self-contained we expose the 
i?-matrix 7?.(A, yu) along with some useful relations among its elements in appendix A. 
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The graded tensor product in 
odromy matrix 

represents the graded Yang-Baxter algebra with the same i?-matrix 



respects comultiphcation. Therefore the mon- 

(45) 



7^(A,/x)(Ti(A)f Ti(M)) = (ri(/x)f T^(A))7^(A,/x) . 



(46) 



We will demonstrate below that Tl{X) generates the Hubbard Hamiltonian. For the 
matrix elements of Tl{X) we introduce the following notation 



V 





Cii 


C12 


D 


fill 


An 


A12 


B 


fi21 


A21 


A22 


B 


D2I 


C21 


C22 


D 



\ 



(47) 



dividing it into four 2x2 submatrices A(A), fi(A), C(A), -D(A). As we shall see in 
sections 5 and 6, this block notation reflects the properties of the monodromy matrix 
under the two su(2) transformations connected with the Hubbard model and under 
combined particle-hole and gauge transformations. 

If a, /5, 7, 5 satisfy (^) and commute among each other and with the diagonal 
elements of Ti(A), then (gT|), (H) and (|U]) imply that 

[tr(rri(A)),tr(rri(/i))] = o . (48) 

Thus tr(r 72(A)) generates a family of mutually commuting operators. Different choices 
of a, /3, 7, 5 correspond to different boundary conditions. Since we did not restrict a, 
/3, 7, 5 to be complex numbers, a dynamical twist is possible. In fact, because of the 
non-local nature of the Jordan- Wigner transformation, the periodic spin model turns 
into the Hubbard model with dynamically twisted boundary conditions. This point was 
recently discussed in detail by Yue and Deguchi |23|. To express the transfer matrix 
t(A) of the spin model introduced in section 2 in terms of fermi operators we note that 
T = Vl^ satisfies (||). Using (0) and (|D we conclude that 



t(A) = ii{v^'rL{\)) 



(49) 
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As we will see in the following, the Hubbard model under periodic boundary conditions 
is obtained with the choice T = cr^ ^ , which leads to 

str(ri(A)) := tr((a^ ® ^')7L(A)) = ti{D) - tr(A) = sti{VLTL{\)) . (50) 

This expression is called the graded trace or super trace of the monodromy matrix. Its 
zeroth order term of the expansion in the spectral parameter is 

str(Ti(0)) = str(ri|(0))str(TL|(0)) = e-'^^^/^^^e-'^'^^/^t/i = e^'^^/^U . (51) 

where N = + is the particle number operator, and U is the shift operator for 
electrons. (|5lD follows from the corresponding result for free fermions which we derive 
in appendix C. U is the product JJ = U^Ui of shift operators for up and down spin 
electrons. We introduce these operators in appendix B, where we give a detailed account 
of shift and momentum operators for fermions on the lattice. If is connected to the total 
momentum 11 by f/ = e'^. 11 assumes its familiar form, 

L-l 

n = 0E^c+5fcr , (52) 
fc=l 

when expressed in terms of Fourier transformed fermi operators, 

gfc. = 4fEe^'^^'Q. . (53) 

1=1 

For brevity we wrote (j) := 27r/L here. Eq. ( ^2]) has to be read with care, since 11 is merely 
defined modulo 27i. Hence we interpret (p2D as defining equation of an equivalence class 
of operators differing from each other only by certain "phase operators". A restriction 
H of H to the fundmental domain of the logarithm is constructed in appendix B. H is a 
polynomial in U. 

The momentum operator H preserves the particle number. Thus 

ln(str(ri(0))) = -i7riV/2 + ifl . (54) 

We will see moreover in section 6 that str(TL(A)) commutes with the particle number 
operator and may therefore conclude that 

[fl,ln(str(ri(A)))] = . (55) 
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This equation implies that r(A) := ln(str(72(A))) is a generating function of transla- 



tional invariant commuting operators. According to the arguments of Liischer |^ these 
operators are local. They are most easily calculated in the language of the spin model, 
since the building blocks of the monodromy matrix are permutation operators of spins, 

To give an example, we present the derivation of the Hamiltonian. Recall from 
section 2 that U = t{0) is the shift operator for spins. If we reintroduce for a moment 
an auxiliary space index a, we obtain 

L-l 

faLiO) = PalULLL-MPLL-1 + E 4+l,(0)P,+l,Palt/ • (56) 

i=i 

The product Ljk{0)Pjk was given in section 2, eq. (0). Using the Jordan- Wigner trans- 
formation, ([27|), (P^), it is expressed in terms of fermi operators as 

Lj + lj iO)Pj + lj = C+ 9 + + + U {uj^ - \){nj^ - \) . (57) 

(0) and (|1|) imply sti{VaLPaiU) = e-'^^I^U. Since moreover [Ki, 4+i,(0)P,+i,] = 0, 
we obtain str(7;L(0)) = e~'''^''^UH. 
H is the Hubbard Hamiltonian 

L 

H = T. (4rC, + lr + C+^lrCjr + f/(n,T - i)(n,i - |)) (58) 

i=i 

under periodic boundary conditions {cl+it '■= Cir). Due to our choice of Jordan- Wigner 
transformation (p7|), (p8| ) we obtained the Hamiltonian for holes here. The sign of 
the hopping term can of course be changed by a particle-hole transformation. Higher 
conserved quantities may be calculated in the same way as the Hamiltonian. We get an 
expansion of the generating function t(A), whose first terms are 

r(A) = -i7riV/2 + ifl + AiJ + C(A2) . (59) 

The (9(A^) term was derived by Shastry 0. The zeroth order terms in ( pQ]) were not 
known before. They are however indispensable for the derivation of the dispersion 
relations of elementary excitations from the eigenvalues of str(7L(A)) [^, |23[. It will 
be interesting to investigate, if we can get both branches of quasiparticle dispersion 
relations P, Q from these eigenvalues. 
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5 Discrete transformations 

A characteristic feature of the Hubbard Hamiltonian on a chain consisting of an even 
number of sites is its invariance under the transformation 

c,T-c,T , c,i-.(-iyc+ , U^-U . (60) 

Since the generators of rotations of the spins are not invariant under (|60|), there exists 
a second su(2) Lie algebra commuting with the Hamiltonian. 

We will show now that not only the Hamiltonian, but the whole transfer matrix 
str(7i(A)) is invariant under (|60D. First note that h{\) —h{\), and thus 

Gh^G-h = Gh^ ■ (61) 

The matrix elements of the free fermion L-matrix (^1]) transform according to 

Al(A) ^e'"*^'^^/V^/:fcl(A)(T^e-'"(^-i)'^^/2 . (62) 

The last two formulae imply 

CkiX) {h ® e'^''^'/') (a' ay)Cki\)ia' ® a^) {h ® e-'^^'-'^'^' , (63) 

and thus by comultiphcation, 

Tl{X) ^ {h ® e'"^'^^/^^ [a' ® ay)TL{X){a' ® a^) . (64) 

Finally str(7L(A)) transforms according to 

str(Ti(A)) ^ -e--^/2 (Dn(A) - e-Mn(A) - A.^iX) + e'^'^D^^iX)) , (65) 

and we can conclude invariance modulo sign of the graded trace for even L 

str(Ti(A)) ^ ±str(r,.(A)) . (66) 

Hence, all higher commuting operators generated by r(A) = ln(str(7^(A))) are invariant 
under the transformation (0). 

We can of course revers the spins in (pOD- Then a slight modification in the trans- 
formation of the monodromy matrix occurs. The factors o"^ ® in (|6^) have to be 
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replaced by a"^ ® I2, and the two factors of I2 ® e''^^""^/^ have to be interchanged. The 
result (^) for the transfer matrix remains the same. Performing up spin and down spin 
transformations in succession the monodromy matrix transforms as 

In this case we find for the transfer matrix 

str(Ti(A)) tr(e-'"^"'Z)(A)) - tr(A(A)) . (68) 
Again invariance is only achieved for an even number of lattice sites. 



6 su(2) symmetries 



A careful discussion of the two su(2) symmetries ||17| , |T8|, |T9|, |2y] connected with the 
Hubbard Hamiltonian was crucial for the complete understanding of the coordinate 
Bethe Ansatz of the model. It was shown in |E5| that the states obtained from coordinate 



Bethe Ansatz are incomplete. They are highest weight states of two su(2) Lie algebras 



p6| , One is the su(2) Lie algebra of rotations, the other one is the ?7-pairing su(2) 
Lie algebra. The generators of the //-pairing su(2) are obtained from the generators 
of rotations under the canonical transformation of the preceding section. They are 
connected with the creation of pairs of particles or holes in the system. 

We show in the following how the generators of the two symmetries commute with 
the monodromy matrix. Our result will be useful for the classification of quasiparticles 
according to their symmetry within the algebraic approach. A discussion analogous to 



the discussion of the spin of spin waves by Faddeev and Takhtajan ||2^ is likely to be 
possible. There are four interactionless states which may serve as reference states for 
an algebraic Bethe Ansatz of the Hubbard model, the empty band, the completely filled 
band and the half-filled band with all spins up or all spins down. Depending on the 
choice of reference state four of the elements of the matrices B{X) and C(A) in ( ^Tj) 
are creation operators, whereas the remaining four are annihilation operators. This fits 
with the fact that there are four different quasiparticles in the system. We think that 
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their identification will eventually become possible by use of the symmetry properties 
presented below. 

The su(2) generators of rotations are given by 

L L L 

j=i j=i j=i 



Recall that we are using the language of holes here (cf. Under a particle hole 

transformation the operators S^, , turn into the operators (\ (, (z used by EBler 
et al. IQ. We will show now that the whole transfer matrix is rotational invariant. To 
this end let us introduce local generators of rotations 

5+ := -c+Cji , Sj- := -c+^Cji , S] := n^^ - n^^ . (70) 

The matrices 

S+:=(T+®(T" , S":=(T"®(T+ , := i((T^®/2 -/2®cT^) (71) 

clearly generate a representation of su(2). They are connected to the inner block A{X) 
of the monodromy matrix Tl{X). Let 

:= + , := -i(S+ - S-) , (72) 

Sj i= Sj ~h Sj 5 S^ i= ) * ('^^) 

Then it is not difficult to see that 

[/:,(A),S" + 5;] =0 , a = x,y,z . (74) 

The verification of this equation may be done as follows. First show by direct calculation 
that [Cj{X), H'^ + Sj'] = 0. Cj{X) has the same block structure as the monodromy matrix, 
(^). Under reversing all spins the blocks A{X), B{X), C{X), D{X) transform as 

/ 







a' 



(75) 



and we may conclude that [Cj{X), S + S'^ ] = 0. The vanishing of the last commutator 
[Cj{X), + 5*1] follows by means of the Jacobi identity. 
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The local equation (j7^) extends to an identity for the monodromy matrix by induc- 
tion, 

[Ti(A),S° + ^"] =0 , (76) 

a = I/, where and are defined as their local analogs. Taking the graded trace 
of this equation yields 

[str(ri(A)),5°] = , a = x,y,z . (77) 

The transfer matrix and thus all higher commuting operators are rotational invariant. 

The transformation properties of the monodromy matrix under the discrete trans- 
formation (^) introduced in the preceding section induce a second su(2) invariance. 
Applying (|60|) to the su(2) generators of rotations, (|69D, we find 



^+ - = E(-l)^'^'44 ' (78) 
i=i 

S- - = E(-l)^-^^c,iC,T , (79) 
i=i 

- ^' = Y^{n,^^n,^-\) = N-L . (80) 
i=i 

This is the //-pairing symmetry. Because of (|80|) , it may be interpreted as non-abelian 
extension of the gauge symmetry. The commutators of the generators of r^-pairing with 
the monodromy matrix follow from (|60D and (|76|). Let 



S+:=a+®cT+ , t-:=(j-®(j- , t"" := \[o' ® h \ h ® o^) . (81) 

These matrices generate a representation of su(2) connected to the block -D(A) of the 
monodromy matrix. Like in case of rotations we define 

:= S+ + S- , := -i(S+ -S") , (82) 

rf := ij^ + T]' , := -i(?7+ - ?7") . (83) 

Using these definitions we obtain for L even 

[ri(A),S" + r7"] = , a = x,y,z , (84) 
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and thus 

[str(ri(A)),r7"]=0 . (85) 
Note that there is no local analog like ( [74| ) of equation (pS]). The ?7-pairing symmetry is 



sensitive to a change of boundary conditions and is in this sense a non-local symmetry. 
Equation (0) may be verified in the following way. Observe that 

J:^{a' 0ay) = {a' 0ay)t^ , (86) 

and 

S±(/2 (g) e'-^^'l'') = (/2 ® e^'"^"'/2)S± . (87) 
Apply the transformation (^) to equation ([76|), and use (p6|), (|87|). Then 

[ri(A),r7±]+Ti(A)S±-(/2®e-^-^)E±r,.(A) = . (88) 

Equation ( |87D remains true, if E='= is replaced by S^. Using this fact, ( pHj) implies 

[TL{X),t^ + V']=0 , (89) 

whereby 

[str(ri(A)),iV] =0 (90) 

for every L, which means that all higher conserved quantities are gauge invariant. This 
fact has been used in the derivation of (|55D . From ( pS]) and (|89D we infer the validity of 
(0) for even L. 

Here is a simple example for the usefulness of the above formulae. (^9]), (^), (|85D 
imply immediately that Ilri~^ = ?7"^(n + 7r), which means that rj^ changes the momentum 



of eigenstates by vr [18 . 



7 Conclusions 

We hope we could convince the reader, that the graded Yang-Baxter algebra ( ^B[ ) is a use- 
ful tool for further investigations of the one-dimensional Hubbard model. Our account 
of basic features of the monodromy matrix should be read in conjunction with the recent 
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preprint of Ramos and Martins |]T6|, who were able to diagonalize str(72(A)) by purely 
algebraic means. Combining both works it should be not too difficult to redrive alge- 
braically all results obtained so far by means of the coordinate Bethe Ansatz. Moreover, 
there is hope to proceed in the calculation of correlation functions and thermodynamical 
properties. 

There have been speculations that there might be a different Yang-Baxter algebra 
embedding of the Hubbard Hamiltonian. Of course we can not rule out this possibility. 
Some arguments in favour of it, however, are disproved by now. In a recent article 



IJ] we were able to show how a rational substructure of the i?-matrix naturally arises 



in the thermodynamic limit. The corresponding sub matrix of the monodromy matrix 
generates the Y(su(2)) representation discovered by Uglov and Korepin [T^. This nicely 
fits with the fact that the ^-matrix of quasiparticle scattering is of rational form. In 
the present article we showed that the monodromy matrix has an appropriate algebraic 
structure. In particular, its graded trace is fully su(2)©su(2) invariant and invariant 
under translations. The analytic properties of i?-matrix and monodromy matrix are less 
usual. We do not yet have a geometrical idea of the spectral curve. Still we succeeded 
in showing the existence of a meromorphic parametrization of the transfer matrix. 
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Appendix A The R-matrix 



The i?-matrix generating the graded Yang-Baxter algebra of the Hubbard model was 
first derived by Olmedilla et al. |Tl|]. It follows from the equations (0), (p3D, ( p^ and 
(1^) and is of the following structure, 

7^(A,/i) = 
/ 
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Pi / 



(A.l) 



The ten Boltzmann weights pj = pj{X,fj,) are 



J-h 



pi = cos(A) cos(yu)e + sin(A) sm{p)e 

P2 = 1 , 

cos(A) cos(/i)e'^~' — sin(A) sin(/i)e'~'* 
cos^(A) — sin^(yu) 

P4 = cos(A) cos{p)e'-~^ + sin(A) sin(/i)e''"' 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
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cos(A) cos(/i)e' ^ — sin(A) sm(fi)e^ 
cos^(A) — sin^(/i) 

2sh(2(/i-/)) 
U{cos^{X) -sm\fi)) 

Pi = Pi- Pb , (Ai 



P8 = pi- P3 , (A.9) 
P9 = sm(A) cos(/x)e'"'' - cos(A) sin(/i)e''~' , (A. 10) 

Pio = sin(A) cos(/i)e''"' - cos(A) sin(/i)e'"'' . (A.ll) 
The parameters A, fi, h and I are connected by equation (P). Note that our definition of 



h and / differs from that in ref. [jTT| and that we performed a shift of ^ in the arguments 
of the functions a{X) and 7(A) occurring there. 



There are the following quadratic relations between the Boltzmann weights, [ p!l[] , 

P1P4 + P9P10 = 1 , (A. 12) 

PiPs + P3P4 = 2 , (A. 13) 

P3P5-P6 = 1 • (A. 14) 

Further identities useful for practical calculations can be found in the recent article 
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Appendix B Momentum operator for fermions on a 
lattice 

Below we present a detailed discussion of the momentum operator on a lattice. 



The shift operator from permutations 

In this subsection we mimic the construction of the shift operator for spin chains. Start 
with spinless fermions, Ci, . . . , c^, on a one dimensional lattice of L sites. Let 

K,,:=l-ict-4)ic,-c,) . (B.l) 
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It is not difficult to see that Kij permutes fermions. There are the obvious identities 

K,, = K+ , K,,=K,, , K,, = l . (B.2) 
Use of the fundamental anti commutators for the fermions yields 

Kijd = Ci + [cf - c^)cj{ci - Cj) = Ci + (-1 + CjC^ - Cjcf){ci - Cj) = CjKij , (B.3) 
and ( |B.2D implies 

Furthermore, by the last four equations, we obtain 

KijKjk = KikKij = KjkKik , j , (B.5) 

and a short calculation similar to the one in eq. (iJ) leads to 

KijKij = l . (B.6) 

Hence, the operators Kij are identified as permutation operators. 

With the aid of permutations of fermions it is of course possible to realize a global 



shift Q. Let 

U:=K^,K23...Kl-il . (B.7) 

Then eq. ( [B.3|) implies 

{Cj+iU j = l,...,L-l 

, if (B.8) 
CiU j=L . 

This means that U is acting as right shift operator on the elementary fermi operators of 
a periodic chain of L sites. Now ( p.7|) implies 

U+ = Kl-il...Ku , (B.9) 

and thus by eq. ( TO , UU+ = U+U = 1. U is invertible, and U'^ = f/+, i.e. U is 
unitary. f/~^ is the left shift operator. With the aid of the shift operator it is possible 
to define the momentum as its formal infinitesimal generator, 

e'" := U . (B.IO) 
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This is the common construction in case of the spin chains. Note however, that ( [B.10|) 
defines 11 only modulo 27c. 

To realize the shift operator for electrons, attach a spin label to all operators above, 
and observe that [f/|,Cj|] = [f/j.,Cj|] = 0. Then U := f/|f/j is the shift operator for 
electrons. 



Diagonalization of the momentum operator — an exercise 

For simplicity, return to the spinless case. Because of the basic property (p.8|) , it follows 
that [U^, Cj] = [U^, Cj'] = 0, j = 1, . . . , L, i.e. is a scalar operator. Since U\0) = |0), 
we obtain = 1. Shifting all fermions once around the lattice does not change the 
state of the system. This simple fact has strong implications. 
Let a E C. Then 

L-l 

(1 - af/) 5] a"f/" = 1 - . (B.ll) 

n=0 

This means that the resolvent of ?7 is a finite sum. With A := l/a we find, 

(A-f/)-i = --^-X:A^-"?7" . (B.12) 

The sum on the right hand side of the latter equation is regular in A and of order A. 
The spectrum of U is therefore given by the equation 

A-^ = l , ^ , Afc = e'2^*^/^ , fc = 0,...,L-l , (B.13) 

and is of course highly degenerate. 

Let Pq, . . . , Pl-1 be the projections on the eigenspaces of U corresponding to the 
eigenvalues Ao,...,Al_i. The are orthogonal, since the eigenstates of a unitary 
operator to different eigenvalues are. PkPi = ^kiPk, and furthermore, Yl,k=Q Pk = 1- 
Hence, the spectral decomposition of the resolvent ( p.l2|) is 



(A - f/r ^ = E tA- ■ (B.14) 

fc=o 
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Conversely, the Pk are determined by the spectral decomposition via = Tesx=x^{{X — 
U)^^). Since d\{X{X^ — l))|A=Afc = L, we obtain 

1 

Pk = YY. e-'^^^'U'' , (B.15) 

^ n=0 

where we used the abbreviation = 2ti/L. From this representation the momentum 
eigenstates arise quite naturally. The particle number operator is of course transla- 
tional invariant, [A^, U] = 0. Hence, the projections Pk preserve the particle number, 
[Pk, N] = 0. We find 

iipjk L 

P,c;|0) = ^Ee-*^"c+|0) . (B.16) 

L n=l 

Since this is true for all j = 1, . . . , L, we see that U is nondegenerate in the one particle 
sector of the Hilbert space. In this sector the subspace corresponding to the eigenvalue 
Afc is spanned by the right hand side of the above equation. Normalization yields the 
eigenvector \k) := c^|0), where is the creator of a one particle momentum eigenstate. 
It is defined in the usual way as the adjoint of the annihilator 

Cfc:=4f Ee'*'"c„ . (B.17) 

Cfc and are fermi operators. They are transformed back into site operators Cj, by 
Fourier inversion of eq. ( |B . 1 7[ ) . Therefore arbitrary operators defined in terms of Cj, Cj 
may be expressed in terms of Ck and . The particle number operator, in particular, is 
form invariant under Fourier transformation, 

N=Y.ctck . (B.18) 

fc=0 

Since the Ck are fermi operators, the states \ki . . . k^) '■= c^^ ■ ■ ■ c^^lO), ki < . . . < kN, 
are orthogonal. ( p.l8|) implies A^|A;i . . . k^) = A^|A;i . . . k^). Counting these states we 
see that they span the A^-particle sector of the lattice Hilbert space. Letting N vary 
from 1, . . . , L we get a basis of the full Hilbert space. 
Now, from the definition ( [B.17| ) of Ck it follows that 



lfc+ = e'^'^ctU . (B.19) 
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Hence we have obtained 

U\ki...kN) = e"t'^^'+-+^''^\h...kN) , (B.20) 

or for the momentum operator, respectively, 

Yi\ki...kN) = <i){ki + ... + kN)\ki...kN) . (B.21) 

This equation is of course defined only modulo 27r. On the other hand it is easily verified 
that the operator J2k=i ^ ^^t^k acts the same way on | fci . . . k^). Since these states form 
a basis, we have achieved the diagonalization of the momentum operator, 

L-1 

Yi = (PY^kctck . (B.22) 
fc=i 

In this form 11 is usually found in the literature. 

Site representation of the momentum operator 

What happens, when we use the Fourier transform to translate back the above diagonal 
form of the momentum operator into the site representation? Inserting the definition 
( PTFQ into eq. yields 

n = ^ E c+ c„ Xi ^e-'^^'"""^' ■ (B.23) 

m,n=l k=l 

Now, for a G C, let 

L—l 1 —i(j>La 

,(a):=Eie-'^'" = iY— =1^ , (B.24) 

k=o ^ ^ 

where the term on the right hand side makes sense only for a ^ LZ. Eqs. ( |B.23D and 
( p:^ imply that 



1 ^ 

n = 7 E 9'{m - n) c+Cn 

m,n=l 



l0(L-l)iV + E .^uiZ:, _ , ■ (B.25) 



m , Ti — 1 



This formula is remarkable. Since 11 = — iln(f/), we can read it as having taken the 
logarithm of the ordered product of transpositions which defines tj . 11 as given by 
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( [B.25| ) is a non-local operator. It is interesting to notice that the second term on the 



right hand side of ( [B.25| ) is gauge equivalent to the 1/ sin(0(m — n)) hopping term of the 
long-range Hubbard model introduced by Gebhard and Ruckenstein |2^ . We have thus 



found a simple interpretation of this Hamiltonian and in turns a simple interpretation 
of the origin of the 1/ sin^ exchange of the Haldane-Shastry spin chain pT| . 



Momentum operator modulo 27r 

To obtain the appropriate definition of 11 mod 27r, observe from the preceding subsection 
that 

-| L L— 1 / 1 iipkm \ 

for k = 0, . . . ,L — 1. With view of the right hand side of this equation the function (pk is 
periodically continued to a saw tooth function on the integers. Since n/0 assumes only 
integer eigenvalues, the definition 

yields the required restriction of 11 modulo 27r. In other words, ( p.l9|) and (p.26| ) imply 
that = e^^. 11 obviously commutes with the Hubbard Hamiltonian, whereas H does 
not. 



Appendix C Free fermions 

In this appendix we show how to obtain the shift operator U for fermions from the 
graded trace of the free fermion monodromy matrix (|33|). For this purpose we have 
to treat XX-chain and free fermion model in parallel. Both models are related by a 
Jordan- Wigner transformation. For brevity let us consider the up-spin case (p7|), and 
let us suppress the spin index here. The XX-chain monodromy matrix is 

TL(A) = /i(A).../i(A) . (C.l) 

It satisfies the Yang-Baxter algebra @, 

f{X-fi){TL{X)®TL{i2)) = {TM®TL{X))f{X-fi) . (C.2) 
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Applying the Jordan- Wigner transformation p7| ) to the XX-chain monodromy matrix, 
we obtain 

Ti(A)=e-'™^'^^/2^i(A) , (C.3) 



where Tl{X) is the free fermion monodromy matrix (p3D with r =|, which satisfies the 
graded analog of ( |C.2| ), 

^.(A-/i)(rz.(A)fr^(/.)) = (r^(^)fr^(A))f,(A-/.) . (c.4) 

The matrix fg{\) is defined as rg{\) ■= W~^f{\)W with W according to (^). ( |C.4|) 
implies that 

[str(ri(A)),str(ri(/i))] = . (C.5) 
We want to calculate the action of str(7i(0)) on fermi operators. To this end let 

Po,:=i(l + aX) , P,fc:=i(l + aX) • (C.6) 
Then (|C.1| ) implies that 

TLiO)=PoiPi2...PL-lL = PoiU . (C.7) 



We use the Jordan- Wigner transformation (pT]) to express Pjj+i in terms of fermi oper- 
ators. Then 

Pjj+i = Kjj+i + 2njnj+i (C.8) 

with the spinless densities Uj = CjCj and the permutation operators Kjj^i according to 
( [BID . Let e± := e^'™^/^ jhen implies that 

str(Ti(0)) = (e+ni + e_(ni - l))U . (C.9) 

Using the Jordan- Wigner transformation (^) and the known action of the permutation 
operators Pjj+i on Pauli matrices, we obtain 

Uck = Ck+i{l-2n^)U , A; = 1,...,L-1 , (C.IO) 
Ucl = Ci{l-2n{)e\U . (C.ll) 
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e+ and e_ are generators of gauge transformations, 

e+Cfc = -icfce+ , e„Cfe = icfcC. . (C.12) 

By use of the four previous formulae, we infer for the annihilators Cfc, 

str(TL(0))cfc = icfc+istr(Ti(0)) , (C.13) 

where k = 1, . . . ,L {cl+i = ci). It is not difficult to see that str(7i(0)) is a unitary 
operator. Therefore the creators satisfy 

str(r^(0))4 = -i4+,str(Ti(0)) , (C.14) 

k = l,...,L. (JfT^) , (^J^) and imply that 



e+str(ri(0)) = «Lf/ , (C.15) 



where is a complex constant, and U is the shift operator (p.Tp for spinless fermions. 
Since Kjj^i\0) = Pjj^i\0) = |0), we conclude from ( |C.9| ) that = — 1 and eventually 
obtain 

str(ri(0)) = -e.U . (C.16) 

The corresponding formula for down-spins is simply obtained by reversing the spins. 
Then (|C.16|) implies (^T]). It may be instructive for the reader to verify ( P.16|) directly 
for small L. 
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